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We report filling-fraction dependent dielectric spectroscopy measurements on the relaxation dy¬ 
namics of the rod-like nematogen 7CB condensed in 13 nm silica nanochannels. In the film-condensed 
regime, a slow interface relaxation dominates the dielectric spectra, whereas from the capillary- 
condensed state up to complete filling an additional, fast relaxation in the core of the channels 
is found. The temperature-dependence of the static capacitance, representative of the averaged, 
collective molecular orientational ordering, indicates a continuous, paranematic-to-nematic (P-N) 
transition, in contrast to the discontinuous bulk behaviour. It is well described by a Landau-de- 
Gennes free energy model for a phase transition in cylindrical conhnement. The large tensile pressure 
of 10 MPa in the capillary-condensed state, resulting from the Young-Laplace pressure at highly 
curved liquid menisci, quantitatively accounts for a downward-shift of the P-N transition and an 
increased molecular mobility in comparison to the unstretched liquid state of the complete filling. 

The strengths of the slow and fast relaxations provide local information on the orientational order: 

The thermotropic behaviour in the core region is bulk-like, i.e. it is characterized by an abrupt onset 
of the nematic order at the P-N transition. By contrast, the interface ordering exhibits a continuous 
evolution at the P-N transition. Thus, the phase behaviour of the entirely filled liquid crystal-silica 
nanocomposite can be quantitatively described by a linear superposition of these distinct nematic 
order contributions. 


I. INTRODUCTION 

Liquid crystals (LCs) spatially confined on the 
nanometer scale exhibit structural and thermodynam¬ 
ical properties which differ often markedly from their 
bulk counterparts. Both the collective orientational 
(isotropic-to-nematic) and the translational (smectic-to- 
liquid or smectic-to-nematic) transitions have turned out 
to be significantly affected by finite size and interfacial 
(solid-liquid or liquid-liquid) interactions introduced by 
confining walls m or the geometrical constraints in 
nanoporous media [H24]. 

For example, the second-order nematic-to-smectic-A 
transition is absent or greatly broadened in rod-like liq¬ 
uid crystals immersed in aerogels. This allowed a de¬ 
tailed study and comparison with theoretical predictions 
of the influence of quenched disorder introduced by ran¬ 
dom spatial confinement on the critical exponents of this 
prominent phase transition m m ng. Similarly, spa¬ 
tial restrictions and interface interactions in mesoporous 
silica and alumina induce pronounced modifications of 
smectic ordering, e.g. often pre-smectic behaviour is ob- 
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servable far above the bulk nematic-to-smectic A tran¬ 
sition. mm It has also been shown experimentally 
[HllillT], in agreement with expectations from theory 
[71 Hg, that there is no ’’true” isotropic-nematic (I-N) 
transition for LCs spatially restricted in at least one di¬ 
rection to a few nanometers. The molecular anchoring 
at the confining walls, quantified by an interfacial field, 
imposes a partial orientational, that results in a partially 
nematic ordering of the confined LCs, even at temper¬ 
atures T far above the bulk I-N transition temperature 
Tjjy. The symmetry breaking and thus entropy change 
does not occur spontaneously, as characteristic of a gen¬ 
uine phase transition. It is enforced over relevant dis¬ 
tances by the interaction with the walls. Thus confine¬ 
ment plays here a similar role as an external magnetic 
field for a spin system [Illig: The strong first order I-N 
transition is replaced by a weak first order or continuous 
paranematic-to-nematic (P-N) transition at a temper¬ 
ature TpN and may also be accompanied by pre transi¬ 
tional phenomena in the molecular orientational distribu¬ 
tion [30j . An understanding of these phenomonenologies 
is of high fundament interest, for it allows to explore the 
validity (and break-down) of basic concepts of condensed 
matter science at the nanoscale [111113 ngug. 

Moreover, LC confined in nanoporous solids provide 
promising composite materials for organic electronics and 
the emerging field of nano photonics, because of the me- 
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FIG. 1: (color online). Schematic side view on distinct meso¬ 
scopic arrangements of a nematic liquid in a nano channel as a 
function of filling fraction / above and below the paranematic- 
to-nematic transition point, Tpm- (a) and (b) refer to the 
film-condensed, (c) and (d) to the capillary-condensed state 
and (e), (f) to the entirely filled channel. The critical filling 
fraction /c separates the film-condensed from the capillary- 
condensed state. The molecules in the interface layer (next 
to the pore walls), at the liquid-vapor interface and inside the 
core region of the channels are marked by different colors and 
labeled as species ”1”, ”2” or ”3”. 


chanical rigidity of the solid host and the large varia¬ 
tion in electrical and optical properties offered by the 
plethora of available LC systems along with the advent 
of tailorable nanoporous media [3TVI34] . 

Referring to the pioneering theoretical work of Sheng, 
Poniewierski, and Sluckin [281 135] as well as experimen¬ 
tal work by Yokovama[36j on LCs in semi-infinite, planar 
confinement, Kutnjak, Kralj, Lahajnar, and Zumer de¬ 
veloped a Landau-de Gennes free energy ansatz [271 El] 
(hereafter denoted as KKLZ model) for the I-N transi¬ 
tion in cylindrical pore geometry. This model describes 
the external orientational field as an effective nematic 
ordering field a. The KKLZ model very successfully 
describes the general features of the P-N transition ob¬ 
served in spatial confinement m mi ESSO], in partic¬ 


ular the T-dependence of the effective order parameter 
averaged over the entire pore geometry. However, it does 
not give any microscopic details on the orientational be¬ 
haviour with regard to spatial inhomogeneities. 

Computer simulations on LCs in thin film and pore ge¬ 
ometry can give here important complementary insights 
mmnnHSH]- These studies indicate pronounced spatial 
heterogeneities, in particular interface-induced molecular 
layering and radial gradients both in the orientational or¬ 
der and reorientational dynamics in cylindrical pore ge¬ 
ometry 1^. Despite recent experimental advancements 
in optical techniques directly probing orientational order 
parameter profiles in the proximity of planar, solid walls 
[52H55] . achieving the spatial (and temporal) resolution 
necessary to rigorously explore such inhomogeneities in 
nanometer-sized capillaries remains still experimentally 
extremely demanding, if not impossible. 

Dielectric spectroscopy on confined LCs, that is a mea¬ 
surement of the frequency-dependent complex electric ca¬ 
pacity C(a;) can make important contributions in that 
respect jSSj. This technique is based on the interaction 
of an external, alternating electric field with the electric 
dipole moments of the molecules. The molecules try to 
arrange their dipoles parallel to the extern field, how¬ 
ever the resulting polarisation is disturbed by thermal 
noise and influenced by the interfacial orientational field. 
The time needed for dipoles to relax and to adjust to 
the alternating electrical field is sensitively depending on 
the local, collective mobility, most prominently the ro¬ 
tational viscosity [SS]. Thus it sensitively affects the di¬ 
electric response, in particular the so-called ” i5-relaxation 
process” which characterizes the movement of the perma¬ 
nent molecular dipole moment about the short geometri¬ 
cal axis of a dielectric positively anisotropic LC-molecule 
[57]. 

A sizeable number of studies on LCs in porous me¬ 
dia evidently document that the rate of dipolar relax¬ 
ations (and thus orientational and translational mobility) 
usually differ significantly between the molecules in the 
pore wall proximity and the ones in the channel center 
[58ff66j . In seminal dielectric studies (complemented by 
dynamic light scattering and calorimetric experiments) 
Aliev et al. [57H5U] could document a slow surface mo¬ 
bility in comparison to the dynamics in the pore centre 
for liquid crystals confined in tortuous and tubular meso- 
pores. Moreover, a significant broadening of the dielectric 
spectra was found and traced by the authors to inhomo¬ 
geneous couplings of the molecules to the pore walls and 
coupling variations among the molecules themselves [70] . 

However, the detailed partitioning, the evolution of the 
relaxation behaviour and the superposition of the dis¬ 
tinct relaxation contributions of LCs as a function of the 
filling of nanoporous media has been elusive so far and 
also controversially discussed gunnia. Here we em¬ 
ploy dielectric relaxation spectroscopy on an archetyp¬ 
ical rod-like liquid crystal (7 CB) confined in an array 
of isolated channels of approximately 7 nm radius in a 
monolithic silica membrane. Characteristic fillings (film 
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condensate, capillary condensate and complete filling) - 
see illustrations in Fig. which we documented in a pre¬ 
vious optical birefringence and light scattering study on 
the identical system [35] , are explored as a function of T. 
This along with the simple parallel capacitor geometry 
(see inset in Fig. allows us to derive detailed insights 
with regard to the heterogeneous mobility and orienta¬ 
tional order of the confined molecular assemblies. Partic¬ 
ularly, we can quantitatively describe the thermotropic 
behaviour with the KKLZ-model described above and 
show that the dielectric response of the organic-inorganic 
nano composite can be described by a linear superposi¬ 
tion of two distinct molecular populations, an interface 
species strongly affected in their molecular mobility and 
orientational order by the silica walls, and the collective 
molecular order in the pore centre. 


II. EXPERIMENTAL 
A. Sample Preparation 

The nematogen liquid crystal 7CB (4-heptyl-4’-cyano 
biphenyl) has been purchased from Merck. The porous 
silicon membranes were prepared by electrochemical an¬ 
odic etching of highly p-doped (100) silicon wafers [73l - 
[75]. The samples were then subjected to thermal oxida¬ 
tion for 12 h at T=800 °C under standard atmosphere. 
The resulting porous silica membranes are permeated 
by parallel channels aligned along the surface normal 
of the membrane m- The average channel radius was 
R = 6.6±0.5 nm (porosity P = 54±2%) according to 
volumetric N 2 -sorption isotherms recorded at T=77 K. 

The entirely filled sample (fraction filling / =1.0) has 
been obtained by melt infiltration, that is capillarity- 
driven filling of the channels [8l|72HZ9|. To prepare the 
partially filled samples (0 < / < 1) we imbibed the 
porous matrix with binary 7CB/cyclohexane solutions of 
selected concentrations. After evaporation of the high- 
vapor-pressure cyclohexane, the matrix remained filled 
by the low-vapor-pressure LC with a filling fraction / 
that depends on the initial solute concentration. The 
filling fraction has been determined by weighting of the 
partially filled matrix and a comparison with the mass 
of the empty matrix. The evaporation of the solvent 
is controlled by time-dependent mass measurements of 
the sample at a temperature of 325 K and followed an 
exponential decay with time. After 1 hour, the cyclo¬ 
hexane had evaporated to 1/100 of its initial concentra¬ 
tion, and thus for the low-/ samples the further mass 
drop was below the detection limits of our balance, which 
was 5 • 10“® g. By a gradual increase of the solute con¬ 
centration, we have prepared three filling fractions, / = 
0.17±0.01, 0.51±0.01 and 0.81±0.01 in addition to the 
completely filled sample, / =1.0. These selected fillings 
cover all known regimes of pore filling, as sketched in 
Fig .[^and documented in a previous optical birefringence 
and light scattering study [55] : the film-condensed regime 
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FIG. 2: (top) Illustration of a dielectric relaxation experiment 
on a monolithic mesoporous membrane traversed by parallel 
nano channels hlled with a liquid. The arrows indicate the 
nano channel’s width and length, respectively, (a) Frequency 
dispersion of the imaginary part of the complex capacitance 
and (b) Cole-Cole plots for four different fractional fillings / 
of 7CB in silica nanochannels {R = 6.6 nm) at T = 295.1 K. 
Symbols correspond to the experimental data points, whereas 
solid lines are hts with a two-relaxation model discussed in the 
text. 


(/ < fc), the capillary-condensed regime {fc < f < 1) 
and the complete filling (/ = 1). Note that according 
to our previous /-dependent optical data sets the criti¬ 
cal filling fraction, separating the film-condensed and the 
capillary-condensed state, is fc » 0.25 [55] . 
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B. Dielectric Relaxation Spectroscopy 

For the dielectric measurements, gold electrodes with 
thicknesses of 50 nm have been sputtered onto the porous 
membrane. All measurements were performed on sam¬ 
ples that had an electrode area of 111 mm^ and a thick¬ 
ness of 340 /im. The geometric capacitance of the samples 
is Co = 2.89 pF. Dielectric spectra have been recorded in 
the frequency range from 0.5 kHz to 15 MHz using the 
Impedance/Gain-Phase Analyzer Solatron-1260A. The 
measurements have been performed at 100 selected tem¬ 
peratures between 292 K and 324 K and the sample tem¬ 
perature was controlled with an accuracy of 0.01 K by a 
LakeShore 340 temperature controller. 


III. RESULTS AND DISCUSSION 

A. Molecular mobility upon channel filling probed 
by dielectric spectroscopy 

In Fig. we display the frequency dependencies of the 
imaginary part of the capacitance C"{v) for four fill¬ 
ing fractions at T =295.1 K. This is far below the P- 
N transition temperature Tpjy = 315 AT, as determined 
in our previous birefringence study [39] . C"{y) is a mea¬ 
sure of the frequency-dependent energy absorption of the 
nanocomposite, when the ^-relaxation occurs. Given the 
absence of any relaxation of the empty silica matrix in 
the frequency regime investigated, it represents a direct 
measure of the molecular relaxation behaviour of the con¬ 
fined LG molecules. Whereas the film-condensed regime 
(/ =0.17) is characterized by a single, broad peak, indi¬ 
cating a slow relaxation with a rather broad distribution 
of relaxation frequencies, another peak, typical of a fast 
relaxation process, evidently emerges on top of its high- 
frequency wing at higher filling fractions, /. It is char¬ 
acterized by a relatively narrow relaxation frequency dis¬ 
tribution and its relaxation strength (height) rises with 
/ up to the entire filling. Thus, the series of dielectric 
dispersions shown in Fig. are quite instructive, since 
they document a partitioning of the pore filling in two 
components with distinct molecular mobility: (i) slowly 
relaxing molecules that are located in a direct contact 
with the pore walls (’’interface species”, see Fig. and 
(ii) fast relaxing molecules that are located in the core 
region of the capillary bridges (a ’’core” component, see 
Fig.0. 

In order to analyse the observed relaxation behaviour 
in more detail, we resort to a representation of the di¬ 
electric data in so-called Cole-Cole diagrams, i.e. we plot 
C"{v) vs. C'{v) with frequency v as the independent 
parameter, see Fig. In such a diagram, a material 
that has a single relaxation time r, as typical of the clas¬ 
sical Debye relaxator, will appear as a semicircle with 
its centre lying on the horizontal (at C" = 0) and the 
peak of the loss factor occurs at l/r. Here, r is a mea¬ 
sure of the dipolar mobility, i. e. it characterises the time 


required by the system to relax by thermally activated 
molecular motions to 1/e times the dipolar ordering in¬ 
duced by the external electrical field. A material with 
a symmetric distribution of relaxation times will be a 
semicircle with its center lying below the horizontal at 
C (y) = 0. An asymmetric distribution of relaxation 
times, e.g. the Havriliak-Negami relaxation, results in an 
asymmetric arc. In that sense, the Cole-Cole representa¬ 
tion allows one to geometrically illustrate and analyse the 
relaxation behaviour of a given system in a quite simple 
manner. [551 [5U] 

As can be seen in Fig. the relaxation behaviour of 
our confined liquid crystalline system is characterized by 
one or two semicircles, depending on the filling fraction. 
Thus, its dielectric relaxation behaviour can be described 
by one or two, so-called Cole-Cole processes: 


^*( 0 ;) 


e*{uj)Co = 

r I I 

°° 1 -k (ia;Ti)i-“i 1 -k (ia;r2)i-“=^ ’ 


( 1 ) 


Here w = 2t:v is the cyclic frequency, e*(uj) is the com¬ 
plex dielectric constant of the composite. Coo = CooCq 
is the high frequency limit capacitance expressed via 
the high frequency permittivity £oo, ACi = AeiCo and 
AC 2 = A£ 2 Co are the capacitance relaxation strengths 
expressed via the dielectric relaxation strengths A£i and 
A £2 of the slow process I and of the fast process H, re¬ 
spectively. The dielectric relaxation strengths quantify 
the amplitude of the different relaxation contributions, 
and thus characterise the extent of dipolar orientational 
polarization typical of each elementary process. The cor¬ 
responding mean relaxation times are ti and T 2 , respec¬ 
tively. The exponent parameter a, can take values be¬ 
tween 0 and 1, and allows one to describe different spec¬ 
tral shapes. When 0 = 0, the Cole-Cole model reduces 
to the Debye model (with a single relaxation frequency). 
When a >0, the relaxation is stretched, i.e. it extends 
over a wider range on a logarithmic frequency scale than 
the Debye process. 

We will see below, that the spectral relaxation shapes 
(and thus a values) change in a characteristic manner as 
a function of temperature and filling fraction for our sys¬ 
tem. However, they are well represented by two symmet¬ 
ric Cole-Cole processes. Deviations of the experimental 
data from the fitting curves at low frequencies can be at¬ 
tributed to ionic DC-conductivity. For this reason, the 
low-frequency region was always excluded from the fitting 
process. A special remark has to be made regarding the 
fitting procedure in the high temperature region (see di¬ 
electric spectra in the electronic supplementary), partic¬ 
ularly just above Tpat, for the fractional fillings, / =0.51, 
0.81 and 1.0. There, the maximum of the imaginary part, 
C''{v) shifts out of the upper limit of the frequency win¬ 
dow {vm >15 MHz) employed in our experiment. Thus 
only a part of the left wing of the relaxation band is ob¬ 
served. This leads to an increased uncertainty of the ex¬ 
tracted fit parameters, especially for the fast relaxation 
process H. Fortunately, one of the parameters affected, 
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namely the capacitance relaxation strength AC' 2 , can be 
determined in an alternative way. Thereby uncertainties 
can be resolved in respect to all other fit parameters for 
this process. In the low-frequency limit {v —>■ 0) Eq. 
gives the static capacitance Cst. It reads: 


Cst = Coo + ACi + AC 2 . (2) 


Since the nanochannels in the silica membrane are 
parallel-aligned to the applied electric field, the entirely 
filled sample forms a simple parallel circuit. Therefore, 
the effective permittivity of the composite is given as 
e = (I - P)esi 02 + CercB, where esiOs = 3.8 and srcB 
are the permittivities of the silica host and 7CB guest, 
respectively. With a Coo =9.5 pF, as determined at 
293 K from a Cole-Cole plot, one obtains Stcb ~ 2.85 
by its extrapolation to high frequencies, i.e. the value 
nearly equal to the dielectric permittivity of 7CB at op¬ 
tical frequencies. This means that the contribution of 
the LC component to the dielectric permittivity at high 
frequencies is due to electronic polarizability only, which 
is weakly temperature dependent. 


Based on refractometric data [ST], we estimated that 
the relative contribution to the changes of the static ca¬ 
pacitance of the composite should not exceed 5%, which 
is within the error margins of our fitting analysis. There¬ 
fore, the assumption of a temperature-independent Coo 
is a good approximation. Moreover, the static dielectric 
constant, Cst, as determined from the Cole-Cole plot by 
its extrapolation to lower frequencies, does not depend on 
other extracted fit parameters. Hence, AC 2 can be alter¬ 
natively evaluated by using the equation]^ The magni¬ 
tudes of AC 2 , representing the extracted fit parameters 
and calculated via equation are depicted in Fig. In 
our fitting procedure the parameter set has been opti¬ 
mised in order to achieve a minimal difference between 
the AC 2 values extracted by this alternative calculation 
procedure. 




FIG. 3: Temperature-dependent capacitance relaxation 
strength of (a) the slow relaxation process I, ACi and (b) 
the fast relaxation process II AC 2 for four selected fraction 
fillings / of 7CB condensed in silica nanochannels. The dash- 
dotted line is the isotropic baselines, ACf°{T) and AC'“°(T) 
(/ =1.0) obtained by the linear extrapolation of AC'i(T) and 
AC 2 {T) from the paranematic region, as discussed in the text. 


Solid lines in Fig. are the best fits obtained by a si¬ 
multaneous analysis of the measured real and imaginary 
parts as outlined above. As a Supplementary the dielec¬ 
tric spectroscopy raw data sets along with fits for several 
selected temperatures and for all four fraction fillings can 
be found. 
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FIG. 4: Temperature-dependent Cole-Cole relaxation distri¬ 
bution parameters of (a) ai, and (b) 02 and the corresponding 
relaxation times (c) n (d), T 2 of the slow and fast relaxation 
process I and II, respectively, as extracted for four character¬ 
istic filling fractions / of silica nanochannels with 7CB. Note 
that the relaxation times are presented in an Arrhenius plot, 
i.e. on a logarithmic scale versus the inverse temperature. 


The relaxation time of the fast process practically co¬ 
incides with that of bulk 7CB [^. It is typical of the 
rotational dynamics around the short molecular axis (5- 
relaxation), the end-over-end tumbling|64[ [68l l82l |83]. 
The T dependences of the extracted fit parameters are 
depicted in Figs. [^and[^ 

In the monolayer regime (/ =0.17) the dielectric relax¬ 
ation is dominated by the slow interface process I. How¬ 


ever, also a weak contribution of the fast relaxation pro¬ 
cess II develops as a function of increasing T - see Fig. IJr. 
This behaviour could originate in nonequivalent positions 
of the nematogen molecules forming the interface layer as 
a function of T. Whereas the filling fraction, / = 0.17 at 
T = 293K, corresponds approximately to one monolayer, 
some molecules are expected to be pushed into the next 
molecular layer at higher Ts, as sketched in Fig. [^, be¬ 
cause of the increase in thermal activity and intermolecu- 
lar spacing. As molecular dynamics studies suggest [44] . 
such molecules are less spatially constrained and thus 
exhibit considerably faster relaxations than those being 
in direct contact with the pore wall (guest-host interac¬ 
tion). Note however that one has to distinguish this fast 
relaxation process in the film-condensed regime from the 
one observed in the capillary-condensed regime or in the 
entirely filled state. Whereas in the first case, the in¬ 
creased mobility is caused by the molecules located at 
the liquid-vapor interface (Fig. type 2 molecules), in 
the second case the fast relaxations are attributable to 
the molecules in the core region of the pore filling with 
negligible pore wall/nematogen interaction (Fig. a type 
3 molecules). This partitioning is corroborated by the 
smooth T variations of all Cole-Cole parameters of the 
relaxation process I and II in the film-condensed regime, 
whereas for higher fractional fillings the relaxation char¬ 
acteristics of the process II are considerably influenced 
by the P-N transition inside the core region, see Figs. 
andjH 

In the capillary condensation regime (/ =0.51 and 
0.81) and for the entire filling (/ =1.0) the Cole-Cole 
parameters, ai and 02 (see Fig.|^) exhibit opposite ten¬ 
dencies in their T variation: The interface relaxation I is 
characterized by a rather broad distribution at room tem¬ 
perature (q:i ~ 0.48-0.50). However, it becomes narrower 
at high Ts and approaches in the paranematic phase a 
value 0.36±0.01. Interestingly, for the fractional Filings, 
/ > 0.51, the corresponding Q:i(r)-dependences are all 
shifted upwards in comparison to the monolayer Filing 
(/ = 0.17). This indicates that by adding molecular lay¬ 
ers one modiFes also the dynamics of the layer in direct 
contact with the silica wall. This observation is intu¬ 
itively understandable and also in good agreement with 
observations on the reorientational dynamics of methanol 
molecules upon silica nanopore Filings |H1|. By contrast, 
the core relaxation H for the entirely Filed sample is 
characterized by a narrow relaxation rate distribution at 
room temperature (a 2 0.06), which broadens in the 

paranematic state, up to a 2 0.11. Moreover, as / 
decreases the 0:2 (r)-curves are systematically shifted up¬ 
wards. One may speculate that this broadening of the re¬ 
laxation time distribution can be attributed to molecules 
located in the menisci regions of the liquid-vapor inter¬ 
faces (see Fig.[^ and d, - type 3 molecules). Their relax¬ 
ation properties should differ from the molecules located 
in the core region, since they encounter more heteroge¬ 
neous surroundings and thus interaction potentials lead¬ 
ing to broader distributions of relaxation times. Also 
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the relative number of these molecules increases with de¬ 
creasing / [39l HSl [86] , in agreement with the increased 
broadening observed for decreasing /. 



FIG. 5: Static capacitance, Cst for four characteristic filling 
fractions / of silica nanochannels with 7CB. The solid line 
is the best fit of the C'st(r)-dependence (/ =1.0) based on 
the KKLZ model. The dash-dotted line is the corresponding 
baseline, obtained by fitting the KKLZ model. 


In order to compare the confined molecular relaxation 
with the documented molecular dynamics of the bulk LC 
system, it is instructive to analyse the T-dependence of 
the relaxation times of the interface and core processes, 
as displayed in Fig. The relaxation times ti{T) and 
T 2 {T) rise strongly with increasing orientational ordering. 
The interface relaxation, however, is much slower than 
the core relaxation process, which can be traced to the at¬ 
tractive silica potential and the expected larger rotational 
and translational viscosities, if one has a continuum pic¬ 
ture in mind [441168) . Both in the nematic and in the 
paranematic phase the fast relaxation in the core region 
exhibit Arrhenius-like behaviour, = Toiexp(Eai/kbT) 
{kb is the Boltzmann constant), typical of a thermally 
activated dipolar reorientation process across an energy 
barrier E^i- For the composites with / =0.51, 0.81 and 
1.0 Arrhenius plots far below Tpjv yields the activation 


energies, Eo 2 , 64.6, 64.9 and 66.2 kJ/mol, respectively, 
for the fast, core i5-relaxation process. These values 
are close to the value of 67.5 kJ/mol reported for bulk 
7CB in Ref |SS|. As expected for the less ordered and 
thus less constrained state, the activation energies are 
smaller in the paranematic phase, namely 39.5, 41.0 and 
43.0 kJ/mol, respectively. Note, however, that in the 
vicinity of Tpjv non-Arrhenius, step-like changes occur 
(see arrows in the Figure). They can be attributed to the 
specific relaxation time behaviour, To 2 , caused by the col¬ 
lective orientational molecular ordering at the P-N tran¬ 
sition. Overall, the behaviour of T 2 {T~^) is analogous to 
the characteristic changes of the static capacitance in the 
vicinity of Tp^v, see Fig.[^ 

The interface relaxation, ti exhibits a smooth variation 
accompanied by a gradual slope change, only. There are 
no sharp features in the vicinity of Tptv- Obviously, the 
thermotropic ordering in the interface layer next to the 
pore wall differs considerably from the one in the core 
region. This conclusion is also corroborated by the be¬ 
haviour of the capacitance relaxation strengths, AC'i(T) 
and AC' 2 (r), as will be discussed in more detail below. 
Comparing the slopes of Arrhenius plots ti{T~^) in the 
composites with different fraction fillings, we again see 
an evident effect of subsequent molecular layers on the 
mobility characteristic of the interface layer. The Ar¬ 
rhenius plot far below Tpjv yields the activation energy, 
Eoi, of 29.8 kJ/mol for the monolayer which rises with / 
reaching a magnitude of 45.4 kJ/mol for the composite 
with / = 1. 


B. Thermotropic nematic order: Landau-de 
Gennes Analysis and Shell/Core Partitioning 

Both relaxation processes contribute to the static 
permittivity with considerably different relaxation 
strengths. As will be outlined in the following this al¬ 
lows us to decompose an effective (average) order param¬ 
eter behaviour in elementary contributions attributable 
to molecular rearrangements in the core and in the shell 
regions, respectively. The orientational order parameter 
Q of a nematic liquid crystal quantifies the collective pre¬ 
ferred molecular orientation in the nematic state. It is 
defined as Q = 5(3cos^ ^ — 1), where 0 is the angle be¬ 
tween the long axis of a 7 CB molecule and the director, 
the mean direction. The brackets mean an averaging 
over the molecules under consideration. In general, Q 
has a local character in a nematic liquid confined in a 
porous medium, since the orientational ordering near the 
pore walls and in the core region may differ |44] . En¬ 
couraged by the observation of the two distinct mobility 
populations, we suggest to approximate the total order 
parameter, Q, which characterizes the molecular order¬ 
ing averaged over the entire pore volume, by a sum of 
two contributions, an interface (or shell), Qs, and a core, 
Qc, contribution: 










Q — Qs “t“ Qc — ^sQs “t“ '^cQc- (3) 

Here the elementary contributions, Qs and Qc, are de¬ 
fined by local order parameters Qs and Qc and weighted 
according to their volume fractions, Ws and Wc with 
Ws + Wc = 1- Note that Qs and Qc are assumed to 
be constant over the distinct pore regions. 

The relation between the static dielectric permittivity 
and the orientational order parameter Q is given by the 
Mayer and Meier equation m- Hence, for a nematic 
molecule having the dipole moment nearly parallel to the 
molecular long axis and with the director aligned parallel 
to the probing field, the capacitance relaxation strength, 
AC can then be expressed in a quite simple way: 

AC = AC^^°+pQ, (4) 

where the first term and the factor of proportionality, p, 
depend on the molecular dipole moment and its orienta¬ 
tion with respect to the long principal axis, the molecu¬ 
lar number density, the internal molecular fields and the 
temperature. Both AC'®® and p are in principle T de¬ 
pendent, but in most cases the corresponding changes 
are weak, particularly if the relative changes in the ab¬ 
solute temperature are small, as is the case here. For 
the bulk nematic LC AC‘®°(T) represents the bare T de¬ 
pendence of the capacitance relaxation strength in the 
isotropic phase, whereas its extrapolation to the nematic 
phase gives an isotropic baseline relative to which an ex¬ 
cess contribution due to orientational ordering (p-term) 
occurs. The Maier-Meier equation can be applied in prin¬ 
ciple to each contribution separately. Thus, the static 
capacitance is given as 

Cst = c'r+p{Qs+Qc) = cr+pQ, ( 5 ) 

where Cf°(r) = Coo + ACf°{T) + ACf°(T) is the 
isotropic baseline of static capacitance, AC'f°{T) and 
AC*’*°(r) are the isotropic baselines of shell and core re¬ 
laxation strengths, respectively. 

Note that in the paranematic state (T > Tpat) Cgt ^ 
Cf°, because of a residual nematic ordering in the inter¬ 
face region. In the presented approach the excess capac¬ 
itance, Cst(T) — Cf°(T) is proportional to the effective 
order parameter, Q, and will be subjected to a further 
analysis. 

The difference, Cst(T) — Cf°iT) (x Q normalized to 
a value at T = 293 K, (5(T)-dependence is depicted in 
Fig. 1^. The effective (averaged) order parameter Q can 
be deconvoluted into elementary contributions character¬ 
izing the molecular orderings in the interface shell {Qs) 
and core {Qc) regions: 

Qs{T)xACi{T)-ACf°{T)- ( 6 ) 

Qc{T) X AC2{T) - A&^°{T). 

where AC'i(T) and AC 2 {T) are the capacitance re¬ 
laxation strengths of the slow and fast processes, re¬ 
spectively. ACf°{T) and ACf°{T) are the isotropic 


baselines. AC‘®°(T) can be easily determined, if one 
takes into account that AC 2 {T) saturates to nearly T 
independent values immediately above the nematic-to- 
paranematic transition point, Tp^- We believe that a lin¬ 
ear extrapolation of this dependence below Tpjv, as it is 
depicted in Fig.[^ (see dash-dot line), properly describes 
the corresponding baseline in the confined nematic phase. 
AC'i(T), on the other hand, evidently exhibits a nonlin¬ 
ear asymptotic behaviour which extents far above Tpjv. 
Nevertheless, the baseline can be determined, if one con¬ 
siders Eq. § ACf°{T) = C'f°{T) - - AC)?°{T), see 

the dash-dot line in Fig. [^a) 

In Fig. we present the T dependences of the static 
capacitance for the composites for all /s investigated. In 
the monolayer regime (/ = 0.17), there is only a grad¬ 
ual increase of Cst upon cooling, indicating a weak in¬ 
crease in orientational ordering. At higher /s pronounced 
changes with a characteristic kink at Tpjv in the Cst{T) 
dependences are observed. Their continuous evolution 
(with a precursor behaviour at high Ts) contrasts with 
the discontinuous behaviour observed at the isotropic-to- 
nematic transition of bulk 7CB m- 

More quantitative insights in this phenomenology can 
be achieved by a comparison with the KKLZ-model 
[271137! for the isotropic-to-nematic transition in cylindri¬ 
cal pore geometries briefly discussed in the introduction. 
According to the KKLZ approach the orientational order¬ 
ing in confinement is characterized by the reduced order 
parameter q = Q/Qo, where Qo = Q(7in) is the bulk 
value of the order parameter taken at the isotropic-to- 
nematic transition temperature. Tin. The dimensionless 
Landau-de Gennes free energy of the nematic LC is then 
represented as 

f = tq^-2q^ + q^ - qa + nq^, (7) 

where t = {T — T*)/ {Tjpf — T*) is the dimensionless tem¬ 
perature [27j . A crucial parameter of the model is the 
difference Tin — T* = AT*, which sets the ’’effective” 
temperature T* and thus calibrates the T-scale of the 
KKLZ-model. Based on the bulk birefringence behaviour 
in the isotropic-to-nematic transition [ID], AT* was de¬ 
termined to 4.9 K. A bilinear coupling between the order 
parameter and the nematic ordering field, i.e. the qa- 
term, is a key feature of the KKLZ theory. It results in 
an upward shift of the effective transition temperature, 
Tppf, and a residual nematic ordering at T > Tpjv, typ¬ 
ical of the paranematic state, as a function of increasing 
strength of the interface ordering field, tr. A cr-value of 
CTc = 0.5 corresponds to a critical point in the a—T phase 
diagram, which separates lines of continuous and discon¬ 
tinuous phase transition evolution The K-term ac¬ 
counts for quenched disorder effects attributable to static 
variations of a, here most prominently attributable to 
pore wall inhomogeneities [HHTj. It rises the free en¬ 
ergy of the low temperature phase and thus provides, in 
contrast to a finite value of a, a downward shift of the 
effective transition temperature. 

The KKLZ model is applicable to entire filled chan- 
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nels only, since the orientational distortions occurring at 
the menisci interfaces for the partially filled state are not 
considered. In addition, our experimental results are af¬ 
fected by depolarization effects at the menisci interfaces. 
Therefore, we restricted our analysis to the entirely filled 
samples (/ = 1). Minimization of Eq. [^with respect to q 
gave the equilibrium value of the scaled order parameter, 
ge, which can be converted to the T-dependence of the 
static capacitance and compared with experiment. By an 
iterative least-square numerical fitting procedure a very 
good agreement between experiment and theory was ob¬ 
tained for ct= 0.76 and k =0.96, see solid curve in Fig. 

It quantitatively accounts for both the paranematic be¬ 
haviour above, and the continuous transition behaviour 
at Tpj^. According to the KKLZ model a is proportional 
to the inverse pore radius, R~^. This allows one to cal¬ 
culate a theoretical, critical channel radius, Rc = Ra/a^ 
m which would separate a continuous (for R < Rc) from 
a discontinuous P-N transition (for R > Rc). It corre¬ 
sponds to 10 nm, a value in agreement with the one ex¬ 
perimentally determined from channel-radius dependent 
birefringence measurements on 7CB |40) . 

Another interesting feature in Fig. is a T-shift of 
the transition temperate, Tpat of AT ~3.5 K, between 
the capillary-condensed fillings (/ =0.51, 0.81) and the 
entirely filled channel (/ =1). One may speculate that 
competing anchoring forces and elastic deformations at 
the liquid-vapor menisci interfaces in the partially filled 
state (as sketched in Fig.[^) could be responsible for this 
stronger preference of the disordered state in the partial 
fillings. There is, however, another, quantitatively ver¬ 
ifiable and arguably more simple explanation, which is 
related to the hydrostatics of the confined liquid: In the 
partially filled, capillary-condensed state, the liquid ex¬ 
periences a tensile pressure, dictated by the concave cur¬ 
vature of the menisci terminating the liquid bridges. This 
negative pressure causes not only subtle deformations of 
the rigid, nanoporous matrix [551150] . it also significantly 
affects density, and thus pressure-dependent phase trans¬ 
formations, most prominently the liquid-solid transition 
[OTHOl] , but also the isotropic-to-nematic transition stud¬ 
ied here [051 IMj ■ According to the Young-Laplace for¬ 
mula applied to the menisci in the capillaries (mean cur¬ 
vature radius = - pore radius = -6.6 nm, surface tension 
of 7 CB=31.7 mN/m [07]) the tensile pressure in the 
liquid bridges amounts to ^-9.6 MPa for /c < / < 1. 
This pressure is completely released upon reaching / = I 
(mean curvature radius = oo). An extrapolation of Tin 
of bulk 7CB as a function of positive, compression pres¬ 
sure, reported in the literature [05], towards this mag¬ 
nitude of negative, tensile pressures yields an expected, 
3.5 K Young-Laplace pressure induced downward shift of 
the transition, in excellent agreement with the AT ob¬ 
served. Therefore, we think that the T-shift is rather an¬ 
other vivid manifestation of the high tensile pressures in 
the capillary-condensed state, than attributable to com¬ 
peting anchoring conditions. This conclusions is further 
corroborated by the observation that the core relaxation 


time T 2 is evidently shorter in the capillary-condensed 
state than in the completely filled state, see Fig. [^ This 
faster molecular reorientation relaxation is also compat¬ 
ible with a negative Young-Laplace pressure. It results 
in a reduced density in the liquid bridges, compared to 
the un-stretched liquid, and thus in a reduced rotational 
viscosity of the confined liquid. 



—,—,—,—,—I—■—I—,—I' 

295 300 305 310 315 320 325 


T[K] 

FIG. 6: The effective (averaged) orientational order parame¬ 
ter, Q vs T (a) and its deconvolution on the elementary local 
order parameters describing molecular orderings in the shell 
{Qa{T)) and core (QdT)) regions (b) - see also the colored 
illustration of the different molecular populations. For conve¬ 
nience, all order parameters are normalized to the value Q at 
293 K. The relation Q = Qs + Qc holds for each T. 


Given the additivity in the Maier-Meier equation for 
the static dielectric relaxations and the resulting orienta¬ 
tional order, we discuss in the following a superposition 
of the orientational order in a shell and in a core con¬ 
tribution under a consideration of a proper treatment 
of the T-dependent baselines of the static capacitance 
contributions (see Methods). The resulting T-dependent 
behaviour of the order parameter averaged over the silica 
nano channels and normalised to its value at T =293 K, 
Q{T) is depicted in Fig. [^. It can be decomposed into 
elementary contributions, calculated by the base-line cor¬ 
rected relaxation strengths, which are characteristic of 
the molecular orderings in the shell {Qs) and core (Qc) 
regions, as displayed in Fig. The T variations of 
the different components in the proximity of Tpiv differ 
considerably, indicating a distinct thermotropic orienta- 
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tional behaviour: Whereas the core ordering exhibits an 
abrupt and almost jump-like increase in nematic order¬ 
ing, reminiscent of the bulk behaviour, the interface or¬ 
dering shows continuous changes with a ’’rounded” kink 
near Tpjv and an asymptotic decrease above this temper¬ 
ature. 

It is understood that this linear decomposition is a 
simplification, for it assumes an entirely decoupled ther¬ 
motropic behaviour between shell and core. There is, 
however, evidence that this is not strictly true: (i) We 
saw changes in the relaxation dynamics of the slow in¬ 
terface process as a function of /. (ii) The core order 
parameter exhibits a small, but systematic decrease with 
T below 305 K. This peculiarity could be an artefact of 
the decomposition, but it also could be real. For ex¬ 
ample, the increase in Qi (and the related stiffening) in 
the interface layer may lead to an increased impact of 
the pore wall irregularities on the thermotropic order in 
the core, and thus to a decrease of Qs. (iii) As one can 
learn from the ACi(T)-dependences, typical of the order 
parameter Qg, the slope of the dependences above Tp^ 
are nearly the same for all /s, see Fig. However, be¬ 
low TpN AC'i(T) remains unchanged for the composite 
with a monolayer (or shell) filling (/ = 0.17), only, and 
obviously exhibits an increasing slope as a function of 
increasing /. This indicates, that the paranematic or¬ 
dering in the host-guest interface evolves independently 
solely for the monolayer filling. By contrast, upon ad¬ 
ditional channel filling the thermotropic ordering in the 
interface layer is enhanced by thermotropic ordering pro¬ 
cesses of the core, presumably via intermolecular interac¬ 
tions. Thus, we believe that the distinct behaviour of the 
relaxation strength and relaxation dynamics (broadening 
and relaxation times) of the two processes well justify a 
qualitative partitioning in a core and a shell behaviour, 
there is, however, still a final interaction between the two 
molecular populations. It is interesting to note that such 
a coupling is also found in Molecular Dynamics simula¬ 
tions by Ji et al. on rod-like LCs in tubular channels [H], 
even though their radially resolved order parameter and 
relaxation time profiles suggest an even stronger mutual 
influence of the shell and the core behaviour as observed 
here. 


IV. CONCLUSIONS 

We reported dielectric studies on the nematic crystal 
7CB confined in parallel-aligned nanochannels in mono¬ 
lithic silica membranes. The measurements have been 
performed on composites characterized by different frac¬ 
tional fillings that cover film condensation at the chan¬ 
nel wall, capillary condensation and the entire filling. 
Whereas for the composite with a LC film a slow inter¬ 
face relaxation dominates the dielectric properties, in all 
other cases the dielectric spectra can be well described by 
two thermally activated Cole-Cole processes with consid¬ 
erably different relaxation rates, i.e. (i) a slow relaxation 


at the interface region in the LC layer next to pore walls 
and (ii) a fast relaxation in the core of the nano channel 
filling. 

Our findings compare well with experimental results 
on the dynamics of ’’conventional” liquid condensates 
[SI ISll I5511103j confined in nanoporous media. In partic¬ 
ular, the observation of distinct core and shell molecular 
mobilities is also in good agreement with experimental 
studies on pore-confined liquid crystalline systems pre¬ 
sented in the past IZIlSSHSi Ei ES], most prominently 
with the work by Aliev et al. discussed in the introduc¬ 
tion. Going beyond this work, we demonstrate here that 
the core ordering is reminiscent of the abrupt, discon¬ 
tinuous bulk order behaviour, whereas the surface or¬ 
dering exhibits a much more continuous evolution with 
a gradual change in slope at the P-N transition point. 
It is asymptotically decreasing above this temperature 
and the superposition of both behaviours results in an ef¬ 
fective isotropic-to-paranematic transition which can be 
well described by a phenomenological Landau-de Gennes 
Ansatz (KKLZ-model). 

These experimental results corroborate the radial par¬ 
titioning and gradual thermotropic orientational order¬ 
ing behaviour found in computer simulations for rod-like 
LCs in cylindrical nanochannels. El m EH [71] Albeit 
the simulation studies indicate a rather smooth transi¬ 
tion in the molecular mobility [H [TT] [72] between the 
interface region and the core, whereas according to our 
experiments this transition is quite sharp, resulting in 
a very distinct monolayer and core behaviour. More¬ 
over, our filling-fraction dependent experiments provide 
information on the thermotropic orientational order as a 
function of mesoscopic arrangement of the nematic liq¬ 
uid in the channels (capillary- and film-condensed states). 
These different fillings states have not been explored in 
simulation studies so far. 

For the future, it would be interesting to what extent 
computer simulations can shed additional light on the 
impact of the tensile pressure in the capillary-condensed 
state in comparison to the influence of possible heteroge¬ 
neous anchoring conditions in this state of the confined 
liquid. In fact, the advent of high brilliant Synchrotron- 
based X-ray sources may allow one to perform already 
in the near future single-channel experiments in order 
to explore spatially resolved anchoring profiles and the 
intimately related translational liquid crystalline order. 
Such structural information would be of high importance 
in order to complement the dynamic and static infor¬ 
mation regarding the orientational order of highly spa¬ 
tially confined nematic liquids unravelled here by dielec¬ 
tric spectroscopy experiments on arrays of nanochannels. 

Supporting Information Available: The frequency dis¬ 
persions of the complex capacitance for seven selected 
temperatures and four fraction fillings of 7CB embedded 
into silica nanochannels {R — 6.6 nm) is available free of 
charge via the Internet as RSC’s Electronic Supplemen¬ 
tary Information (ESI). 
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